home *** CD-ROM | disk | FTP | other *** search
/ Enter 2004 January / enter-2004-01.iso / files / maxima-5.9.0.exe / {app} / share / maxima / 5.9.0 / src / numerical / slatec / dyairy.lisp < prev    next >
Encoding:
Text File  |  2003-02-09  |  30.4 KB  |  565 lines
  1. ;;; Compiled by f2cl version 2.0 beta 2002-05-06
  2. ;;; 
  3. ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
  4. ;;;           (:coerce-assigns :as-needed) (:array-type ':simple-array)
  5. ;;;           (:array-slicing nil) (:declare-common nil)
  6. ;;;           (:float-format double-float))
  7.  
  8. (in-package "SLATEC")
  9.  
  10.  
  11. (let ((n1 20)
  12.       (n2 19)
  13.       (n3 14)
  14.       (m1 18)
  15.       (m2 17)
  16.       (m3 12)
  17.       (n1d 21)
  18.       (n2d 20)
  19.       (n3d 19)
  20.       (n4d 14)
  21.       (m1d 19)
  22.       (m2d 18)
  23.       (m3d 17)
  24.       (m4d 12)
  25.       (fpi12 1.30899693899575)
  26.       (spi12 1.83259571459405)
  27.       (con1 0.666666666666667)
  28.       (con2 7.74148278841779)
  29.       (con3 0.36476610549035604)
  30.       (bk1 (make-array 20 :element-type 'double-float))
  31.       (bk2 (make-array 20 :element-type 'double-float))
  32.       (bk3 (make-array 20 :element-type 'double-float))
  33.       (bk4 (make-array 14 :element-type 'double-float))
  34.       (bjp (make-array 19 :element-type 'double-float))
  35.       (bjn (make-array 19 :element-type 'double-float))
  36.       (aa (make-array 14 :element-type 'double-float))
  37.       (bb (make-array 14 :element-type 'double-float))
  38.       (dbk1 (make-array 21 :element-type 'double-float))
  39.       (dbk2 (make-array 20 :element-type 'double-float))
  40.       (dbk3 (make-array 20 :element-type 'double-float))
  41.       (dbk4 (make-array 14 :element-type 'double-float))
  42.       (dbjp (make-array 19 :element-type 'double-float))
  43.       (dbjn (make-array 19 :element-type 'double-float))
  44.       (daa (make-array 14 :element-type 'double-float))
  45.       (dbb (make-array 14 :element-type 'double-float)))
  46.   (declare (type (simple-array double-float (21)) dbk1)
  47.            (type (simple-array double-float (19)) dbjn dbjp bjn bjp)
  48.            (type (simple-array double-float (14)) dbb daa dbk4 bb aa bk4)
  49.            (type (simple-array double-float (20)) dbk3 dbk2 bk3 bk2 bk1)
  50.            (type double-float con3 con2 con1 spi12 fpi12)
  51.            (type f2cl-lib:integer4 m4d m3d m2d m1d n4d n3d n2d n1d m3 m2 m1 n3
  52.             n2 n1))
  53.   (f2cl-lib:fset (f2cl-lib:fref bk1 (1) ((1 20))) 2.43202846447449)
  54.   (f2cl-lib:fset (f2cl-lib:fref bk1 (2) ((1 20))) 2.57132009754685)
  55.   (f2cl-lib:fset (f2cl-lib:fref bk1 (3) ((1 20))) 1.02802341258616)
  56.   (f2cl-lib:fset (f2cl-lib:fref bk1 (4) ((1 20))) 0.34195817820587204)
  57.   (f2cl-lib:fset (f2cl-lib:fref bk1 (5) ((1 20))) 0.0841978629889284)
  58.   (f2cl-lib:fset (f2cl-lib:fref bk1 (6) ((1 20))) 0.0193877282587962)
  59.   (f2cl-lib:fset (f2cl-lib:fref bk1 (7) ((1 20))) 0.00392687837130335)
  60.   (f2cl-lib:fset (f2cl-lib:fref bk1 (8) ((1 20))) 6.83302689948043e-4)
  61.   (f2cl-lib:fset (f2cl-lib:fref bk1 (9) ((1 20))) 1.14611403991141e-4)
  62.   (f2cl-lib:fset (f2cl-lib:fref bk1 (10) ((1 20))) 1.7419513833708603e-5)
  63.   (f2cl-lib:fset (f2cl-lib:fref bk1 (11) ((1 20))) 2.41223620956355e-6)
  64.   (f2cl-lib:fset (f2cl-lib:fref bk1 (12) ((1 20))) 3.2452559198327297e-7)
  65.   (f2cl-lib:fset (f2cl-lib:fref bk1 (13) ((1 20))) 4.03509798540183e-8)
  66.   (f2cl-lib:fset (f2cl-lib:fref bk1 (14) ((1 20))) 4.70875059642296e-9)
  67.   (f2cl-lib:fset (f2cl-lib:fref bk1 (15) ((1 20))) 5.35367432585889e-10)
  68.   (f2cl-lib:fset (f2cl-lib:fref bk1 (16) ((1 20))) 5.70606721846334e-11)
  69.   (f2cl-lib:fset (f2cl-lib:fref bk1 (17) ((1 20))) 5.80526363709933e-12)
  70.   (f2cl-lib:fset (f2cl-lib:fref bk1 (18) ((1 20))) 5.763389886163881e-13)
  71.   (f2cl-lib:fset (f2cl-lib:fref bk1 (19) ((1 20))) 5.42103834518071e-14)
  72.   (f2cl-lib:fset (f2cl-lib:fref bk1 (20) ((1 20))) 4.91857330301677e-15)
  73.   (f2cl-lib:fset (f2cl-lib:fref bk2 (1) ((1 20))) 0.574830555784088)
  74.   (f2cl-lib:fset (f2cl-lib:fref bk2 (2) ((1 20))) -0.006916486483768911)
  75.   (f2cl-lib:fset (f2cl-lib:fref bk2 (3) ((1 20))) 0.00197460263052093)
  76.   (f2cl-lib:fset (f2cl-lib:fref bk2 (4) ((1 20))) -5.24043043868823e-4)
  77.   (f2cl-lib:fset (f2cl-lib:fref bk2 (5) ((1 20))) 1.2296514723966104e-4)
  78.   (f2cl-lib:fset (f2cl-lib:fref bk2 (6) ((1 20))) -2.2705951446217298e-5)
  79.   (f2cl-lib:fset (f2cl-lib:fref bk2 (7) ((1 20))) 2.23575555008526e-6)
  80.   (f2cl-lib:fset (f2cl-lib:fref bk2 (8) ((1 20))) 4.1517495502389895e-7)
  81.   (f2cl-lib:fset (f2cl-lib:fref bk2 (9) ((1 20))) -2.8498575219823097e-7)
  82.   (f2cl-lib:fset (f2cl-lib:fref bk2 (10) ((1 20))) 8.501871747754352e-8)
  83.   (f2cl-lib:fset (f2cl-lib:fref bk2 (11) ((1 20))) -1.7040082689132602e-8)
  84.   (f2cl-lib:fset (f2cl-lib:fref bk2 (12) ((1 20))) 2.2547974674688903e-9)
  85.   (f2cl-lib:fset (f2cl-lib:fref bk2 (13) ((1 20))) -1.0952416657744299e-10)
  86.   (f2cl-lib:fset (f2cl-lib:fref bk2 (14) ((1 20))) -3.4106384509971094e-11)
  87.   (f2cl-lib:fset (f2cl-lib:fref bk2 (15) ((1 20))) 1.11262893886662e-11)
  88.   (f2cl-lib:fset (f2cl-lib:fref bk2 (16) ((1 20))) -1.75542944241734e-12)
  89.   (f2cl-lib:fset (f2cl-lib:fref bk2 (17) ((1 20))) 1.36298600401767e-13)
  90.   (f2cl-lib:fset (f2cl-lib:fref bk2 (18) ((1 20))) 8.763421057556641e-15)
  91.   (f2cl-lib:fset (f2cl-lib:fref bk2 (19) ((1 20))) -4.64063099157041e-15)
  92.   (f2cl-lib:fset (f2cl-lib:fref bk2 (20) ((1 20))) 7.7877275873296e-16)
  93.   (f2cl-lib:fset (f2cl-lib:fref bk3 (1) ((1 20))) 0.566777053506912)
  94.   (f2cl-lib:fset (f2cl-lib:fref bk3 (2) ((1 20))) 0.0026367282834957903)
  95.   (f2cl-lib:fset (f2cl-lib:fref bk3 (3) ((1 20))) 5.1230335147313e-5)
  96.   (f2cl-lib:fset (f2cl-lib:fref bk3 (4) ((1 20))) 2.10229231564492e-6)
  97.   (f2cl-lib:fset (f2cl-lib:fref bk3 (5) ((1 20))) 1.4221709511389e-7)
  98.   (f2cl-lib:fset (f2cl-lib:fref bk3 (6) ((1 20))) 1.2853429589126403e-8)
  99.   (f2cl-lib:fset (f2cl-lib:fref bk3 (7) ((1 20))) 7.28556219407507e-10)
  100.   (f2cl-lib:fset (f2cl-lib:fref bk3 (8) ((1 20))) -3.45236157301011e-10)
  101.   (f2cl-lib:fset (f2cl-lib:fref bk3 (9) ((1 20))) -2.11919115912724e-10)
  102.   (f2cl-lib:fset (f2cl-lib:fref bk3 (10) ((1 20))) -6.56803892922376e-11)
  103.   (f2cl-lib:fset (f2cl-lib:fref bk3 (11) ((1 20))) -8.14873160315074e-12)
  104.   (f2cl-lib:fset (f2cl-lib:fref bk3 (12) ((1 20))) 3.0317784563218303e-12)
  105.   (f2cl-lib:fset (f2cl-lib:fref bk3 (13) ((1 20))) 1.73447220554115e-12)
  106.   (f2cl-lib:fset (f2cl-lib:fref bk3 (14) ((1 20))) 1.67935548701554e-13)
  107.   (f2cl-lib:fset (f2cl-lib:fref bk3 (15) ((1 20))) -1.49622868806719e-13)
  108.   (f2cl-lib:fset (f2cl-lib:fref bk3 (16) ((1 20))) -5.154704589534069e-14)
  109.   (f2cl-lib:fset (f2cl-lib:fref bk3 (17) ((1 20))) 8.757418418578301e-15)
  110.   (f2cl-lib:fset (f2cl-lib:fref bk3 (18) ((1 20))) 7.9673555352572e-15)
  111.   (f2cl-lib:fset (f2cl-lib:fref bk3 (19) ((1 20))) -1.29566137861742e-16)
  112.   (f2cl-lib:fset (f2cl-lib:fref bk3 (20) ((1 20))) -1.1187879441752e-15)
  113.   (f2cl-lib:fset (f2cl-lib:fref bk4 (1) ((1 14))) 0.48544438670511403)
  114.   (f2cl-lib:fset (f2cl-lib:fref bk4 (2) ((1 14))) -0.0030852508840846305)
  115.   (f2cl-lib:fset (f2cl-lib:fref bk4 (3) ((1 14))) 6.987484048379281e-5)
  116.   (f2cl-lib:fset (f2cl-lib:fref bk4 (4) ((1 14))) -2.8275723417976795e-6)
  117.   (f2cl-lib:fset (f2cl-lib:fref bk4 (5) ((1 14))) 1.59553313064138e-7)
  118.   (f2cl-lib:fset (f2cl-lib:fref bk4 (6) ((1 14))) -1.12980692144601e-8)
  119.   (f2cl-lib:fset (f2cl-lib:fref bk4 (7) ((1 14))) 9.476715154987541e-10)
  120.   (f2cl-lib:fset (f2cl-lib:fref bk4 (8) ((1 14))) -9.083017360264231e-11)
  121.   (f2cl-lib:fset (f2cl-lib:fref bk4 (9) ((1 14))) 9.707762064507239e-12)
  122.   (f2cl-lib:fset (f2cl-lib:fref bk4 (10) ((1 14))) -1.1368752725457401e-12)
  123.   (f2cl-lib:fset (f2cl-lib:fref bk4 (11) ((1 14))) 1.4398291753341502e-13)
  124.   (f2cl-lib:fset (f2cl-lib:fref bk4 (12) ((1 14))) -1.9521101955881498e-14)
  125.   (f2cl-lib:fset (f2cl-lib:fref bk4 (13) ((1 14))) 2.81056379909357e-15)
  126.   (f2cl-lib:fset (f2cl-lib:fref bk4 (14) ((1 14))) -4.26916444775176e-16)
  127.   (f2cl-lib:fset (f2cl-lib:fref bjp (1) ((1 19))) 0.134918611457638)
  128.   (f2cl-lib:fset (f2cl-lib:fref bjp (2) ((1 19))) -0.31931458820581304)
  129.   (f2cl-lib:fset (f2cl-lib:fref bjp (3) ((1 19))) 0.0522061946276114)
  130.   (f2cl-lib:fset (f2cl-lib:fref bjp (4) ((1 19))) 0.0528869112170312)
  131.   (f2cl-lib:fset (f2cl-lib:fref bjp (5) ((1 19))) -0.0085810075607735)
  132.   (f2cl-lib:fset (f2cl-lib:fref bjp (6) ((1 19))) -0.00299211002025555)
  133.   (f2cl-lib:fset (f2cl-lib:fref bjp (7) ((1 19))) 4.21126741969759e-4)
  134.   (f2cl-lib:fset (f2cl-lib:fref bjp (8) ((1 19))) 8.739318303692731e-5)
  135.   (f2cl-lib:fset (f2cl-lib:fref bjp (9) ((1 19))) -1.0674916347753302e-5)
  136.   (f2cl-lib:fset (f2cl-lib:fref bjp (10) ((1 19))) -1.5657509725934898e-6)
  137.   (f2cl-lib:fset (f2cl-lib:fref bjp (11) ((1 19))) 1.6805115198399898e-7)
  138.   (f2cl-lib:fset (f2cl-lib:fref bjp (12) ((1 19))) 1.89901103638691e-8)
  139.   (f2cl-lib:fset (f2cl-lib:fref bjp (13) ((1 19))) -1.81374004961922e-9)
  140.   (f2cl-lib:fset (f2cl-lib:fref bjp (14) ((1 19))) -1.66339134593739e-10)
  141.   (f2cl-lib:fset (f2cl-lib:fref bjp (15) ((1 19))) 1.4295633578080996e-11)
  142.   (f2cl-lib:fset (f2cl-lib:fref bjp (16) ((1 19))) 1.10179811626595e-12)
  143.   (f2cl-lib:fset (f2cl-lib:fref bjp (17) ((1 19))) -8.60187724192263e-14)
  144.   (f2cl-lib:fset (f2cl-lib:fref bjp (18) ((1 19))) -5.71248177285064e-15)
  145.   (f2cl-lib:fset (f2cl-lib:fref bjp (19) ((1 19))) 4.084145528538029e-16)
  146.   (f2cl-lib:fset (f2cl-lib:fref bjn (1) ((1 19))) 0.06590416735256971)
  147.   (f2cl-lib:fset (f2cl-lib:fref bjn (2) ((1 19))) -0.424905910566004)
  148.   (f2cl-lib:fset (f2cl-lib:fref bjn (3) ((1 19))) 0.28720974519583)
  149.   (f2cl-lib:fset (f2cl-lib:fref bjn (4) ((1 19))) 0.12978777109960601)
  150.   (f2cl-lib:fset (f2cl-lib:fref bjn (5) ((1 19))) -0.0456354317590358)
  151.   (f2cl-lib:fset (f2cl-lib:fref bjn (6) ((1 19))) -0.010263017598254001)
  152.   (f2cl-lib:fset (f2cl-lib:fref bjn (7) ((1 19))) 0.00250704671521101)
  153.   (f2cl-lib:fset (f2cl-lib:fref bjn (8) ((1 19))) 3.78127183743483e-4)
  154.   (f2cl-lib:fset (f2cl-lib:fref bjn (9) ((1 19))) -7.112875832840839e-5)
  155.   (f2cl-lib:fset (f2cl-lib:fref bjn (10) ((1 19))) -8.08651210688923e-6)
  156.   (f2cl-lib:fset (f2cl-lib:fref bjn (11) ((1 19))) 1.2387953127328499e-6)
  157.   (f2cl-lib:fset (f2cl-lib:fref bjn (12) ((1 19))) 1.13096815867279e-7)
  158.   (f2cl-lib:fset (f2cl-lib:fref bjn (13) ((1 19))) -1.4623428317631e-8)
  159.   (f2cl-lib:fset (f2cl-lib:fref bjn (14) ((1 19))) -1.11576315688077e-9)
  160.   (f2cl-lib:fset (f2cl-lib:fref bjn (15) ((1 19))) 1.24846618243897e-10)
  161.   (f2cl-lib:fset (f2cl-lib:fref bjn (16) ((1 19))) 8.18334132555274e-12)
  162.   (f2cl-lib:fset (f2cl-lib:fref bjn (17) ((1 19))) -8.071748770484839e-13)
  163.   (f2cl-lib:fset (f2cl-lib:fref bjn (18) ((1 19))) -4.6377861876642495e-14)
  164.   (f2cl-lib:fset (f2cl-lib:fref bjn (19) ((1 19))) 4.09043399081631e-15)
  165.   (f2cl-lib:fset (f2cl-lib:fref aa (1) ((1 14))) -0.278593552803079)
  166.   (f2cl-lib:fset (f2cl-lib:fref aa (2) ((1 14))) 0.00352915691882584)
  167.   (f2cl-lib:fset (f2cl-lib:fref aa (3) ((1 14))) 2.31149677384994e-5)
  168.   (f2cl-lib:fset (f2cl-lib:fref aa (4) ((1 14))) -4.7131784226356e-6)
  169.   (f2cl-lib:fset (f2cl-lib:fref aa (5) ((1 14))) 1.12415907931333e-7)
  170.   (f2cl-lib:fset (f2cl-lib:fref aa (6) ((1 14))) 2.0010030118433902e-8)
  171.   (f2cl-lib:fset (f2cl-lib:fref aa (7) ((1 14))) -2.6094807530219305e-9)
  172.   (f2cl-lib:fset (f2cl-lib:fref aa (8) ((1 14))) 3.55098136101216e-11)
  173.   (f2cl-lib:fset (f2cl-lib:fref aa (9) ((1 14))) 3.50849978423875e-11)
  174.   (f2cl-lib:fset (f2cl-lib:fref aa (10) ((1 14))) -5.830071879542019e-12)
  175.   (f2cl-lib:fset (f2cl-lib:fref aa (11) ((1 14))) 2.04644828753326e-13)
  176.   (f2cl-lib:fset (f2cl-lib:fref aa (12) ((1 14))) 1.1052917947674198e-13)
  177.   (f2cl-lib:fset (f2cl-lib:fref aa (13) ((1 14))) -2.87724778038775e-14)
  178.   (f2cl-lib:fset (f2cl-lib:fref aa (14) ((1 14))) 2.8820511100993906e-15)
  179.   (f2cl-lib:fset (f2cl-lib:fref bb (1) ((1 14))) -0.490275424742791)
  180.   (f2cl-lib:fset (f2cl-lib:fref bb (2) ((1 14))) -0.00157647277946204)
  181.   (f2cl-lib:fset (f2cl-lib:fref bb (3) ((1 14))) 9.66195963140306e-5)
  182.   (f2cl-lib:fset (f2cl-lib:fref bb (4) ((1 14))) -1.35916080268815e-7)
  183.   (f2cl-lib:fset (f2cl-lib:fref bb (5) ((1 14))) -2.98157342654859e-7)
  184.   (f2cl-lib:fset (f2cl-lib:fref bb (6) ((1 14))) 1.8682476755997898e-8)
  185.   (f2cl-lib:fset (f2cl-lib:fref bb (7) ((1 14))) 1.0368573766714102e-9)
  186.   (f2cl-lib:fset (f2cl-lib:fref bb (8) ((1 14))) -3.28660818434328e-10)
  187.   (f2cl-lib:fset (f2cl-lib:fref bb (9) ((1 14))) 2.5709141063278e-11)
  188.   (f2cl-lib:fset (f2cl-lib:fref bb (10) ((1 14))) 2.3235765530067695e-12)
  189.   (f2cl-lib:fset (f2cl-lib:fref bb (11) ((1 14))) -9.575232790482551e-13)
  190.   (f2cl-lib:fset (f2cl-lib:fref bb (12) ((1 14))) 1.2034082804971902e-13)
  191.   (f2cl-lib:fset (f2cl-lib:fref bb (13) ((1 14))) 2.9090771677071503e-15)
  192.   (f2cl-lib:fset (f2cl-lib:fref bb (14) ((1 14))) -4.5565645458014903e-15)
  193.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (1) ((1 21))) 2.95926143981893)
  194.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (2) ((1 21))) 3.86774568440103)
  195.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (3) ((1 21))) 1.80441072356289)
  196.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (4) ((1 21))) 0.5780707641253281)
  197.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (5) ((1 21))) 0.16301146817470802)
  198.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (6) ((1 21))) 0.0392044409961855)
  199.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (7) ((1 21))) 0.007909642104338121)
  200.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (8) ((1 21))) 0.0015064086316733801)
  201.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (9) ((1 21))) 2.56651976920042e-4)
  202.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (10) ((1 21))) 3.93826605867715e-5)
  203.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (11) ((1 21))) 5.810977714638179e-6)
  204.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (12) ((1 21))) 7.86881233754659e-7)
  205.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (13) ((1 21))) 9.93272957325739e-8)
  206.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (14) ((1 21))) 1.2142420557510702e-8)
  207.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (15) ((1 21))) 1.38528332697707e-9)
  208.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (16) ((1 21))) 1.50190067586758e-10)
  209.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (17) ((1 21))) 1.58271945457594e-11)
  210.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (18) ((1 21))) 1.57531847699042e-12)
  211.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (19) ((1 21))) 1.50774055398181e-13)
  212.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (20) ((1 21))) 1.40594335806564e-14)
  213.   (f2cl-lib:fset (f2cl-lib:fref dbk1 (21) ((1 21))) 1.2494269877721802e-15)
  214.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (1) ((1 20))) 0.5497568094324711)
  215.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (2) ((1 20))) 0.00913556983276901)
  216.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (3) ((1 20))) -0.00253635048605507)
  217.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (4) ((1 20))) 6.604237953420541e-4)
  218.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (5) ((1 20))) -1.5521724313541602e-4)
  219.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (6) ((1 20))) 3.00090325448633e-5)
  220.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (7) ((1 20))) -3.7645433946734796e-6)
  221.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (8) ((1 20))) -1.3329133161161597e-7)
  222.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (9) ((1 20))) 2.4258737104901298e-7)
  223.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (10) ((1 20))) -8.07861075240228e-8)
  224.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (11) ((1 20))) 1.7109281886119299e-8)
  225.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (12) ((1 20))) -2.4108735757059901e-9)
  226.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (13) ((1 20))) 1.53910848162371e-10)
  227.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (14) ((1 20))) 2.5646537319062995e-11)
  228.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (15) ((1 20))) -9.88581911653212e-12)
  229.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (16) ((1 20))) 1.6087798641263099e-12)
  230.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (17) ((1 20))) -1.2095252474173898e-13)
  231.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (18) ((1 20))) -1.0697827841082e-14)
  232.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (19) ((1 20))) 5.024785570675611e-15)
  233.   (f2cl-lib:fset (f2cl-lib:fref dbk2 (20) ((1 20))) -8.689861309358862e-16)
  234.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (1) ((1 20))) 0.560598509354302)
  235.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (2) ((1 20))) -0.00364870013248135)
  236.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (3) ((1 20))) -5.98147152307417e-5)
  237.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (4) ((1 20))) -2.3361159525362496e-6)
  238.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (5) ((1 20))) -1.64571516521436e-7)
  239.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (6) ((1 20))) -2.06333012920569e-8)
  240.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (7) ((1 20))) -4.2774543157311e-9)
  241.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (8) ((1 20))) -1.08494137799276e-9)
  242.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (9) ((1 20))) -2.3720718887276304e-10)
  243.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (10) ((1 20))) -2.22132920864966e-11)
  244.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (11) ((1 20))) 1.0723800803213801e-11)
  245.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (12) ((1 20))) 5.719548452458081e-12)
  246.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (13) ((1 20))) 7.511027377778351e-13)
  247.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (14) ((1 20))) -3.8191236948379303e-13)
  248.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (15) ((1 20))) -1.7587005711925702e-13)
  249.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (16) ((1 20))) 6.69641694419084e-15)
  250.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (17) ((1 20))) 2.26866724792055e-14)
  251.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (18) ((1 20))) 2.69898141356743e-15)
  252.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (19) ((1 20))) -2.6713361239735905e-15)
  253.   (f2cl-lib:fset (f2cl-lib:fref dbk3 (20) ((1 20))) -6.54121403165269e-16)
  254.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (1) ((1 14))) 0.49307299918803604)
  255.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (2) ((1 14))) 0.00438335419803815)
  256.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (3) ((1 14))) -8.37413882246205e-5)
  257.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (4) ((1 14))) 3.2026881048463196e-6)
  258.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (5) ((1 14))) -1.7566197954827e-7)
  259.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (6) ((1 14))) 1.22269906524508e-8)
  260.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (7) ((1 14))) -1.01381314366052e-9)
  261.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (8) ((1 14))) 9.636397842374749e-11)
  262.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (9) ((1 14))) -1.02344993379648e-11)
  263.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (10) ((1 14))) 1.19264576554355e-12)
  264.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (11) ((1 14))) -1.5044389910328698e-13)
  265.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (12) ((1 14))) 2.0329905237934903e-14)
  266.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (13) ((1 14))) -2.9189065200829206e-15)
  267.   (f2cl-lib:fset (f2cl-lib:fref dbk4 (14) ((1 14))) 4.4232208197547496e-16)
  268.   (f2cl-lib:fset (f2cl-lib:fref dbjp (1) ((1 19))) 0.113140872390745)
  269.   (f2cl-lib:fset (f2cl-lib:fref dbjp (2) ((1 19))) -0.208301511416328)
  270.   (f2cl-lib:fset (f2cl-lib:fref dbjp (3) ((1 19))) 0.0169396341953138)
  271.   (f2cl-lib:fset (f2cl-lib:fref dbjp (4) ((1 19))) 0.0290895212478621)
  272.   (f2cl-lib:fset (f2cl-lib:fref dbjp (5) ((1 19))) -0.0034146713131154903)
  273.   (f2cl-lib:fset (f2cl-lib:fref dbjp (6) ((1 19))) -0.00146455339197417)
  274.   (f2cl-lib:fset (f2cl-lib:fref dbjp (7) ((1 19))) 1.6331327289851702e-4)
  275.   (f2cl-lib:fset (f2cl-lib:fref dbjp (8) ((1 19))) 3.9114532892216197e-5)
  276.   (f2cl-lib:fset (f2cl-lib:fref dbjp (9) ((1 19))) -3.96757190808119e-6)
  277.   (f2cl-lib:fset (f2cl-lib:fref dbjp (10) ((1 19))) -6.51846913772395e-7)
  278.   (f2cl-lib:fset (f2cl-lib:fref dbjp (11) ((1 19))) 5.9870749526928e-8)
  279.   (f2cl-lib:fset (f2cl-lib:fref dbjp (12) ((1 19))) 7.441086545365491e-9)
  280.   (f2cl-lib:fset (f2cl-lib:fref dbjp (13) ((1 19))) -6.212410565226321e-10)
  281.   (f2cl-lib:fset (f2cl-lib:fref dbjp (14) ((1 19))) -6.18768017313526e-11)
  282.   (f2cl-lib:fset (f2cl-lib:fref dbjp (15) ((1 19))) 4.7232348475232394e-12)
  283.   (f2cl-lib:fset (f2cl-lib:fref dbjp (16) ((1 19))) 3.9165245980253194e-13)
  284.   (f2cl-lib:fset (f2cl-lib:fref dbjp (17) ((1 19))) -2.7498593784522596e-14)
  285.   (f2cl-lib:fset (f2cl-lib:fref dbjp (18) ((1 19))) -1.9503649776275003e-15)
  286.   (f2cl-lib:fset (f2cl-lib:fref dbjp (19) ((1 19))) 1.26669643809444e-16)
  287.   (f2cl-lib:fset (f2cl-lib:fref dbjn (1) ((1 19))) -0.018809126006885)
  288.   (f2cl-lib:fset (f2cl-lib:fref dbjn (2) ((1 19))) -0.14779818082614002)
  289.   (f2cl-lib:fset (f2cl-lib:fref dbjn (3) ((1 19))) 0.546075900433171)
  290.   (f2cl-lib:fset (f2cl-lib:fref dbjn (4) ((1 19))) 0.15214693266311602)
  291.   (f2cl-lib:fset (f2cl-lib:fref dbjn (5) ((1 19))) -0.0958260412266886)
  292.   (f2cl-lib:fset (f2cl-lib:fref dbjn (6) ((1 19))) -0.016310273169613)
  293.   (f2cl-lib:fset (f2cl-lib:fref dbjn (7) ((1 19))) 0.00575364806680105)
  294.   (f2cl-lib:fset (f2cl-lib:fref dbjn (8) ((1 19))) 7.12145408252655e-4)
  295.   (f2cl-lib:fset (f2cl-lib:fref dbjn (9) ((1 19))) -1.7545211684672402e-4)
  296.   (f2cl-lib:fset (f2cl-lib:fref dbjn (10) ((1 19))) -1.71063171685128e-5)
  297.   (f2cl-lib:fset (f2cl-lib:fref dbjn (11) ((1 19))) 3.2443558063168e-6)
  298.   (f2cl-lib:fset (f2cl-lib:fref dbjn (12) ((1 19))) 2.61190663932884e-7)
  299.   (f2cl-lib:fset (f2cl-lib:fref dbjn (13) ((1 19))) -4.03026865912779e-8)
  300.   (f2cl-lib:fset (f2cl-lib:fref dbjn (14) ((1 19))) -2.7643516585389505e-9)
  301.   (f2cl-lib:fset (f2cl-lib:fref dbjn (15) ((1 19))) 3.5968792906231206e-10)
  302.   (f2cl-lib:fset (f2cl-lib:fref dbjn (16) ((1 19))) 2.14953308456051e-11)
  303.   (f2cl-lib:fset (f2cl-lib:fref dbjn (17) ((1 19))) -2.4184931190390102e-12)
  304.   (f2cl-lib:fset (f2cl-lib:fref dbjn (18) ((1 19))) -1.2806800492075102e-13)
  305.   (f2cl-lib:fset (f2cl-lib:fref dbjn (19) ((1 19))) 1.2693983440177303e-14)
  306.   (f2cl-lib:fset (f2cl-lib:fref daa (1) ((1 14))) 0.277571356944231)
  307.   (f2cl-lib:fset (f2cl-lib:fref daa (2) ((1 14))) -0.0044421283341992005)
  308.   (f2cl-lib:fset (f2cl-lib:fref daa (3) ((1 14))) 8.42328522190089e-5)
  309.   (f2cl-lib:fset (f2cl-lib:fref daa (4) ((1 14))) 2.5804031841870995e-6)
  310.   (f2cl-lib:fset (f2cl-lib:fref daa (5) ((1 14))) -3.4238972021762104e-7)
  311.   (f2cl-lib:fset (f2cl-lib:fref daa (6) ((1 14))) 6.24286894709776e-9)
  312.   (f2cl-lib:fset (f2cl-lib:fref daa (7) ((1 14))) 2.3637783684457703e-9)
  313.   (f2cl-lib:fset (f2cl-lib:fref daa (8) ((1 14))) -3.16991042656673e-10)
  314.   (f2cl-lib:fset (f2cl-lib:fref daa (9) ((1 14))) 4.4099569165819097e-12)
  315.   (f2cl-lib:fset (f2cl-lib:fref daa (10) ((1 14))) 5.186742210935751e-12)
  316.   (f2cl-lib:fset (f2cl-lib:fref daa (11) ((1 14))) -9.64874015137022e-13)
  317.   (f2cl-lib:fset (f2cl-lib:fref daa (12) ((1 14))) 4.9019057660870996e-14)
  318.   (f2cl-lib:fset (f2cl-lib:fref daa (13) ((1 14))) 1.77253430678112e-14)
  319.   (f2cl-lib:fset (f2cl-lib:fref daa (14) ((1 14))) -5.55950610442662e-15)
  320.   (f2cl-lib:fset (f2cl-lib:fref dbb (1) ((1 14))) 0.49162732110460106)
  321.   (f2cl-lib:fset (f2cl-lib:fref dbb (2) ((1 14))) 0.00311164930427489)
  322.   (f2cl-lib:fset (f2cl-lib:fref dbb (3) ((1 14))) 8.23140762854081e-5)
  323.   (f2cl-lib:fset (f2cl-lib:fref dbb (4) ((1 14))) -4.6176977617214193e-6)
  324.   (f2cl-lib:fset (f2cl-lib:fref dbb (5) ((1 14))) -6.13158880534626e-8)
  325.   (f2cl-lib:fset (f2cl-lib:fref dbb (6) ((1 14))) 2.8729580465652e-8)
  326.   (f2cl-lib:fset (f2cl-lib:fref dbb (7) ((1 14))) -1.8195971537211703e-9)
  327.   (f2cl-lib:fset (f2cl-lib:fref dbb (8) ((1 14))) -1.44752826642035e-10)
  328.   (f2cl-lib:fset (f2cl-lib:fref dbb (9) ((1 14))) 4.5372404342042194e-11)
  329.   (f2cl-lib:fset (f2cl-lib:fref dbb (10) ((1 14))) -3.99655065847223e-12)
  330.   (f2cl-lib:fset (f2cl-lib:fref dbb (11) ((1 14))) -3.24089119830323e-13)
  331.   (f2cl-lib:fset (f2cl-lib:fref dbb (12) ((1 14))) 1.6209895256874102e-13)
  332.   (f2cl-lib:fset (f2cl-lib:fref dbb (13) ((1 14))) -2.4076524797405702e-14)
  333.   (f2cl-lib:fset (f2cl-lib:fref dbb (14) ((1 14))) 1.6938481128449098e-16)
  334.   (defun dyairy (x rx c bi dbi)
  335.     (declare (type double-float dbi bi c rx x))
  336.     (prog ((ax 0.0) (cv 0.0) (d1 0.0) (d2 0.0) (ex 0.0) (e1 0.0) (e2 0.0)
  337.            (f1 0.0) (f2 0.0) (rtrx 0.0) (s1 0.0) (s2 0.0) (t_ 0.0) (tc 0.0)
  338.            (temp1 0.0) (temp2 0.0) (tt 0.0) (i 0) (j 0))
  339.       (declare (type f2cl-lib:integer4 j i)
  340.                (type double-float tt temp2 temp1 tc t_ s2 s1 rtrx f2 f1 e2 e1
  341.                 ex d2 d1 cv ax))
  342.       (setf ax (coerce (abs x) 'double-float))
  343.       (setf rx (f2cl-lib:fsqrt ax))
  344.       (setf c (* con1 ax rx))
  345.       (if (< x 0.0) (go label120))
  346.       (if (> c 8.0) (go label60))
  347.       (if (> x 2.5) (go label30))
  348.       (setf t_ (* (- (+ x x) 2.5) 0.4))
  349.       (setf tt (+ t_ t_))
  350.       (setf j n1)
  351.       (setf f1 (f2cl-lib:fref bk1 (j) ((1 20))))
  352.       (setf f2 0.0)
  353.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  354.                     ((> i m1) nil)
  355.         (tagbody
  356.           (setf j (f2cl-lib:int-sub j 1))
  357.           (setf temp1 f1)
  358.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref bk1 (j) ((1 20)))))
  359.           (setf f2 temp1)
  360.          label10))
  361.       (setf bi (+ (- (* t_ f1) f2) (f2cl-lib:fref bk1 (1) ((1 20)))))
  362.       (setf j n1d)
  363.       (setf f1 (f2cl-lib:fref dbk1 (j) ((1 21))))
  364.       (setf f2 0.0)
  365.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  366.                     ((> i m1d) nil)
  367.         (tagbody
  368.           (setf j (f2cl-lib:int-sub j 1))
  369.           (setf temp1 f1)
  370.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dbk1 (j) ((1 21)))))
  371.           (setf f2 temp1)
  372.          label20))
  373.       (setf dbi (+ (- (* t_ f1) f2) (f2cl-lib:fref dbk1 (1) ((1 21)))))
  374.       (go end_label)
  375.      label30
  376.       (setf rtrx (f2cl-lib:fsqrt rx))
  377.       (setf t_ (* (- (+ x x) con2) con3))
  378.       (setf tt (+ t_ t_))
  379.       (setf j n1)
  380.       (setf f1 (f2cl-lib:fref bk2 (j) ((1 20))))
  381.       (setf f2 0.0)
  382.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  383.                     ((> i m1) nil)
  384.         (tagbody
  385.           (setf j (f2cl-lib:int-sub j 1))
  386.           (setf temp1 f1)
  387.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref bk2 (j) ((1 20)))))
  388.           (setf f2 temp1)
  389.          label40))
  390.       (setf bi (/ (+ (- (* t_ f1) f2) (f2cl-lib:fref bk2 (1) ((1 20)))) rtrx))
  391.       (setf ex (exp c))
  392.       (setf bi (* bi ex))
  393.       (setf j n2d)
  394.       (setf f1 (f2cl-lib:fref dbk2 (j) ((1 20))))
  395.       (setf f2 0.0)
  396.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  397.                     ((> i m2d) nil)
  398.         (tagbody
  399.           (setf j (f2cl-lib:int-sub j 1))
  400.           (setf temp1 f1)
  401.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dbk2 (j) ((1 20)))))
  402.           (setf f2 temp1)
  403.          label50))
  404.       (setf dbi
  405.               (* (+ (- (* t_ f1) f2) (f2cl-lib:fref dbk2 (1) ((1 20)))) rtrx))
  406.       (setf dbi (* dbi ex))
  407.       (go end_label)
  408.      label60
  409.       (setf rtrx (f2cl-lib:fsqrt rx))
  410.       (setf t_ (- (/ 16.0 c) 1.0))
  411.       (setf tt (+ t_ t_))
  412.       (setf j n1)
  413.       (setf f1 (f2cl-lib:fref bk3 (j) ((1 20))))
  414.       (setf f2 0.0)
  415.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  416.                     ((> i m1) nil)
  417.         (tagbody
  418.           (setf j (f2cl-lib:int-sub j 1))
  419.           (setf temp1 f1)
  420.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref bk3 (j) ((1 20)))))
  421.           (setf f2 temp1)
  422.          label70))
  423.       (setf s1 (+ (- (* t_ f1) f2) (f2cl-lib:fref bk3 (1) ((1 20)))))
  424.       (setf j n2d)
  425.       (setf f1 (f2cl-lib:fref dbk3 (j) ((1 20))))
  426.       (setf f2 0.0)
  427.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  428.                     ((> i m2d) nil)
  429.         (tagbody
  430.           (setf j (f2cl-lib:int-sub j 1))
  431.           (setf temp1 f1)
  432.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dbk3 (j) ((1 20)))))
  433.           (setf f2 temp1)
  434.          label80))
  435.       (setf d1 (+ (- (* t_ f1) f2) (f2cl-lib:fref dbk3 (1) ((1 20)))))
  436.       (setf tc (+ c c))
  437.       (setf ex (exp c))
  438.       (if (> tc 35.0) (go label110))
  439.       (setf t_ (- (/ 10.0 c) 1.0))
  440.       (setf tt (+ t_ t_))
  441.       (setf j n3)
  442.       (setf f1 (f2cl-lib:fref bk4 (j) ((1 14))))
  443.       (setf f2 0.0)
  444.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  445.                     ((> i m3) nil)
  446.         (tagbody
  447.           (setf j (f2cl-lib:int-sub j 1))
  448.           (setf temp1 f1)
  449.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref bk4 (j) ((1 14)))))
  450.           (setf f2 temp1)
  451.          label90))
  452.       (setf s2 (+ (- (* t_ f1) f2) (f2cl-lib:fref bk4 (1) ((1 14)))))
  453.       (setf bi (/ (+ s1 (* (exp (- tc)) s2)) rtrx))
  454.       (setf bi (* bi ex))
  455.       (setf j n4d)
  456.       (setf f1 (f2cl-lib:fref dbk4 (j) ((1 14))))
  457.       (setf f2 0.0)
  458.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  459.                     ((> i m4d) nil)
  460.         (tagbody
  461.           (setf j (f2cl-lib:int-sub j 1))
  462.           (setf temp1 f1)
  463.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dbk4 (j) ((1 14)))))
  464.           (setf f2 temp1)
  465.          label100))
  466.       (setf d2 (+ (- (* t_ f1) f2) (f2cl-lib:fref dbk4 (1) ((1 14)))))
  467.       (setf dbi (* rtrx (+ d1 (* (exp (- tc)) d2))))
  468.       (setf dbi (* dbi ex))
  469.       (go end_label)
  470.      label110
  471.       (setf bi (/ (* ex s1) rtrx))
  472.       (setf dbi (* ex rtrx d1))
  473.       (go end_label)
  474.      label120
  475.       (if (> c 5.0) (go label150))
  476.       (setf t_ (- (* 0.4 c) 1.0))
  477.       (setf tt (+ t_ t_))
  478.       (setf j n2)
  479.       (setf f1 (f2cl-lib:fref bjp (j) ((1 19))))
  480.       (setf e1 (f2cl-lib:fref bjn (j) ((1 19))))
  481.       (setf f2 0.0)
  482.       (setf e2 0.0)
  483.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  484.                     ((> i m2) nil)
  485.         (tagbody
  486.           (setf j (f2cl-lib:int-sub j 1))
  487.           (setf temp1 f1)
  488.           (setf temp2 e1)
  489.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref bjp (j) ((1 19)))))
  490.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref bjn (j) ((1 19)))))
  491.           (setf f2 temp1)
  492.           (setf e2 temp2)
  493.          label130))
  494.       (setf bi
  495.               (- (+ (- (* t_ e1) e2) (f2cl-lib:fref bjn (1) ((1 19))))
  496.                  (* ax (+ (- (* t_ f1) f2) (f2cl-lib:fref bjp (1) ((1 19)))))))
  497.       (setf j n3d)
  498.       (setf f1 (f2cl-lib:fref dbjp (j) ((1 19))))
  499.       (setf e1 (f2cl-lib:fref dbjn (j) ((1 19))))
  500.       (setf f2 0.0)
  501.       (setf e2 0.0)
  502.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  503.                     ((> i m3d) nil)
  504.         (tagbody
  505.           (setf j (f2cl-lib:int-sub j 1))
  506.           (setf temp1 f1)
  507.           (setf temp2 e1)
  508.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dbjp (j) ((1 19)))))
  509.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref dbjn (j) ((1 19)))))
  510.           (setf f2 temp1)
  511.           (setf e2 temp2)
  512.          label140))
  513.       (setf dbi
  514.               (+ (* x x (+ (- (* t_ f1) f2) (f2cl-lib:fref dbjp (1) ((1 19)))))
  515.                  (+ (- (* t_ e1) e2) (f2cl-lib:fref dbjn (1) ((1 19))))))
  516.       (go end_label)
  517.      label150
  518.       (setf rtrx (f2cl-lib:fsqrt rx))
  519.       (setf t_ (- (/ 10.0 c) 1.0))
  520.       (setf tt (+ t_ t_))
  521.       (setf j n3)
  522.       (setf f1 (f2cl-lib:fref aa (j) ((1 14))))
  523.       (setf e1 (f2cl-lib:fref bb (j) ((1 14))))
  524.       (setf f2 0.0)
  525.       (setf e2 0.0)
  526.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  527.                     ((> i m3) nil)
  528.         (tagbody
  529.           (setf j (f2cl-lib:int-sub j 1))
  530.           (setf temp1 f1)
  531.           (setf temp2 e1)
  532.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref aa (j) ((1 14)))))
  533.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref bb (j) ((1 14)))))
  534.           (setf f2 temp1)
  535.           (setf e2 temp2)
  536.          label160))
  537.       (setf temp1 (+ (- (* t_ f1) f2) (f2cl-lib:fref aa (1) ((1 14)))))
  538.       (setf temp2 (+ (- (* t_ e1) e2) (f2cl-lib:fref bb (1) ((1 14)))))
  539.       (setf cv (- c fpi12))
  540.       (setf bi (/ (+ (* temp1 (cos cv)) (* temp2 (sin cv))) rtrx))
  541.       (setf j n4d)
  542.       (setf f1 (f2cl-lib:fref daa (j) ((1 14))))
  543.       (setf e1 (f2cl-lib:fref dbb (j) ((1 14))))
  544.       (setf f2 0.0)
  545.       (setf e2 0.0)
  546.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  547.                     ((> i m4d) nil)
  548.         (tagbody
  549.           (setf j (f2cl-lib:int-sub j 1))
  550.           (setf temp1 f1)
  551.           (setf temp2 e1)
  552.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref daa (j) ((1 14)))))
  553.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref dbb (j) ((1 14)))))
  554.           (setf f2 temp1)
  555.           (setf e2 temp2)
  556.          label170))
  557.       (setf temp1 (+ (- (* t_ f1) f2) (f2cl-lib:fref daa (1) ((1 14)))))
  558.       (setf temp2 (+ (- (* t_ e1) e2) (f2cl-lib:fref dbb (1) ((1 14)))))
  559.       (setf cv (- c spi12))
  560.       (setf dbi (* (- (* temp1 (cos cv)) (* temp2 (sin cv))) rtrx))
  561.       (go end_label)
  562.      end_label
  563.       (return (values nil rx c bi dbi)))))
  564.  
  565.